Shear banding and flow-concentration coupling in colloidal glasses 
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We report experiments on hard sphere colloidal glasses that show a type of shear banding hitherto 
unobserved in soft glasses. We present a scenario that relates this to an instability due to shear- 
concentration coupling, a mechanism previously thought unimportant in these materials. Below a 
characteristic shear rate 7 C we observe increasingly non-linear and localized velocity profiles. We 
attribute this to very slight concentration gradients in the unstable flow regime. A simple model 
accounts for both the observed increase of y c with concentration, and the fluctuations in the flow. 
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Shear banding is widespread in the flow of disordered 
materials, including complex fluids [l|; pastes, gels and 
emulsions Q ; granular matter [|[ , soils and rocks Q and 
metallic glasses Q. Understanding this phenomenon is 
thus crucial in various fields of science and engineering. 

Constitutive models of shear banding exist, e.g. in 
wormlike micellar fluids [l|, where a flow curve of stress 
versus strain rate cr( A f) with slope rid < causes instabil- 
ity and separation into two bands with distinct flow rates 
(7i > 72) > [l|; these can also have distinct concentra- 
tions [6]. Many systems with a yield stress (a y > 0) 
also show coexistence of distinct bands with one 

band now being solid (71 =0). In some cases this stems 
from a similar mechanical instability < for 7 < 72) 
f9l-fl2} due to positive feedback between flow and struc- 
tural breakup. However the banding seen in other exper- 
iments [TH O, [H, [13] and simulations [H, Eil , cannot 
be explained in those terms, particularly for purely repul- 
sive interactions [2(| ■ In some of these systems, banding 
may be attributed to cooperativity between local plastic 
events, characterized by a cooperativity length £. How- 
ever, theory [2l| and simulations [22J show that £ grows 
when the flow slows down, which contrasts with the rate- 
independence seen in experiments [lfjj . 

In this Letter we show by experiment and theory that 
concentrated hard-sphere (HS) colloids, one of the sim- 
plest yield-stress fluids and a model for soft glasses gener- 
ally, can exhibit a type of shear banding that does not fit 
into any of the above categories. Instead, we propose a 
scenario where banding is caused by shear-concentration 
coupling (SCC). Though well known as a generic mecha- 
nism for flow instability [23j , this has not previously been 
explored as a shear banding mechanism in glasses. This 
is perhaps because the concentration changes involved 
can be extremely small, as we show below; hence they 
are not directly detectable in experiments. Crucially, the 
effects on flow of very small concentration gradients are 
vastly amplified by the presence of a yield stress. 

It is well known that for nonuniform stress, particle mi- 
gration [24[ is driven by gradients in 7 and the nonequi- 
librium particle pressure, 11(7) [25], [26|. The resulting 
concentration inhomogeneity causes departures from flow 



profiles for the homogeneous system. But this can not 
explain such departures under uniform stress, which re- 
quires intrinsic instability (23[ (see also [Hj|). Recent the- 
ories |27j show that Newtonian and linearly viscoelastic 



materials can also exhibit such instability via SCC. How- 
ever, this arises at 45° to the flow direction, and is unre- 
lated to the results in [23[ on SCC-induced instability in 
nonlinear fluids. In applying the latter to glasses, our two 
key ingredients are (non-linear) dilatancy -the tendency 
of jammed systems to expand under flow, (dH/dj > 0) 
and flow nonlinearity. As both are ubiquitous in glassy 
materials, our results are likely to have wide relevance. 

We used sterically stabilised polymethylmethacrylate 
particles (radii a — 138 nm and 150 nm, polydispersity 
~ 15%) suspended in a decalin-tetralin mixture (viscos- 
ity r/ s — 2.3 mPas) and seeded with ~ 0.5% fluorescent 
colloids (a — 652 nm) of the same kind. Different vol- 
ume fractions <f> were prepared by diluting samples cen- 
trifuged to a sediment with cf> — (f> m ~ 0.67; we report 
data in terms of the reduced concentration $ = 4>/<f> m . 
The glass transition was found to lie at 3> 9 ~ 0.86 [28j. 

Rheology was measured in an AR2000 rheometer in 
cone/transparent plate geometry (cone angle 1°, radius 
20 mm), coupled to a confocal microscope [2{| [3(| to 
measure the velocity v(z) across the gap, < z < z g , at 
various z g . For 7 > 0.1 s _1 we measured v(z) at each 
z for a time > 3/7; for 7 < 0.1 s _1 , v(z) was measured 
from rapid 3D scans [30( over a time ~ 2/7. The reported 
v(z) are constant over the measurement time. To prevent 
slip [29| . walls were coated with a disordered monolayer 
of the tracers; results for a cone and/or plate roughened 
to ~ 10 fj,m were very similar. We took data at fixed 7, 
stepping down 7 after 300 s preshear at 7 ~ 10 — 30 s _1 or 
stepping up from 7 = (waiting for > 2/7 before acquir- 
ing v(z) in either case). The results showed no system- 
atic differences between these protocols; stress-controlled 
measurements also showed no significant changes. 

For $ > $ g , the bulk rheology is as previously re- 
ported, Fig. [Ua), with flow curves of Herschcl-Bulklcy 
(HB) form: a - a y oc 7", with n ~ 0.4 - 0.5. The strong 
increase of cr v with <&, due to the vanishing of free volume 



as $ 1 [32|, is consistent with a y (§) ~ cr (l - $) 
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FIG. 1: (a) Flow curves with HB fits (lines) for $ = 0.875 
and <& = 0.933. (b) Yield stress versus <&. Line: a y = ao(l — 
$)~ 3 with a = 0Mk B T/a 3 . (c) Velocity profiles v(z) for 
$ = 0.933 at various 7. (d) Evolution of v(z) after startup 
shear of 7 = 0.01 s" 1 at to ($ = 0.94); data for t-to = 3740 s 
show error bars. Data in (a),(c) and (d) are for a = 138 nm. 



with p ~ 3 and cr ~ 0.01fc B T/a 3 , Fig. [TJb). 

Until now, HB and similar monotonic flow curves for 
glasses have not been linked to non-transient shear band- 
ing. However, the underlying velocity profiles, shown in 
Fig. QJc) for z g = 170 /im and $ = 0.933, exhibit a 
marked change when we decrease the imposed shear rate 
7 = J A f(z)dz/ z g (here 7(2) = d z v). At large 7, v(z) 
is linear, but for 7 = 0.2 s _1 v(z) becomes highly non- 
linear, with an enhanced rate near the plate and a pro- 
gressive reduction towards the cone. For even smaller 
rates, 7 < 0.05 s" 1 , the nonlinearity grows and 7(2) de- 
creases continuously from a value 3> 7 near the plate to 
7(2) ~ for larger z. The width of the fluidized band 
appears to saturate for low 7 at ~ 80a (not shown); we 
find no evidence for a minimum strain rate in this band 
[3l| . The behavior for z g = 60 /im and 90 /im is essen- 
tially the same. Such continuously varying flow profiles 
strongly contrast with the distinct solid and fluid bands 
in thixotropic yield stress fluids H,[I3|. HS glasses, which 
show only very weak aging of quiescent properties [32T ]. 
are thus distinct from such systems. Note from Figs.[]7c) 
and [2jb) that v(z) has no unique 'symmetry': the flu- 
idized band may appear near either the plate or the cone. 
This rules out sedimentation or specific wall rheology [l6| 
as explanations [33|. The growth of the fluidized band 
with 7 also contrasts with the rate dependence (or lack 
thereof) of the cooperativity length found in [l6ll2lT|. 

Next, we discuss the concentration dependence of the 
observed behavior, Fig. [2j For both <£> = 0.895 (just 
within the glass, Fig. HJa)) and for a much higher con- 
centration $ = 0.948, Fig.^b)), we again observe a tran- 
sition to a nonlinear velocity profile as 7 is lowered, but 
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FIG. 2: Velocity profiles at different volume fractions: (a) 
$ = 0.895 (a = 138 nm), (b) $ = 0.948 (a = 150 nm). 



the shear rate at which this occurs is respectively much 
smaller and higher than for $ = 0.933, Fig. [ljc). On 
decreasing 7, we define the critical shear rate, j c , to be 
that point at which the maximum deviation from linear- 
ity of the normalized velocity profile, \v(z) — 7z|/v C onc, 
first exceeds 0.1. Results for the critical Peclet number, 
Pe c = 7 c tb (where tb — 6irrj s a 3 /ksT is the Brownian 
time) are shown as a function of $ in Fig. |3l 

To begin to interpret our observations, we first show 
that the observed velocity profiles can be reconciled with 
HB behavior, simply by postulating a small concentra- 
tion variation 5§(z) across the gap. Writing the HB form 
as yr = [{a / a v {<S>)) - l] 1 /™ with o y {>S>) = a Q (l - $_)- 3 
as before, we can then calculate v(z, a) = 7(<E> + 
SQ(z'),a)dz' for a given mean concentration <t and a 
choice of <5<£>(z). In Fig. BJa) we do this for a uniform 
gradient d5<S>(z)/dz = \6$\/z g with |<5$|/$_= 0.002, at 
various values of the reduced stress o/o y (<&). When a 
approaches o~ y (&), v(z) changes from weakly to strongly 
nonlinear, reflecting the progressive localization of shear 
within regions of the sample with the lowest yield stress 
tj y ($(z)), i.e. with the lowest $(2). These results strik- 
ingly resemble the experimental data in Fig. [lib), al- 
though <5<£> is too small to be directly measured [34j; dif- 
ferent symmetries of v(z) in other experiments can also 
be explained by corresponding changes in 8$(z). Note 
that the mean shear rate 7 = v(z g , a, <$>{z))/z g differs 
from 7 (cr, $), but the effective flow curves (7(7, &(z)) de- 
viate only slightly from the uniform a(j, $), see Fig.QJb). 

While concentration gradients can thus account for the 
results, their origin and the enhanced shear localization 
with increasing 3> remain to be explained. We now show 
that both are explicable via the SCC instability scenario 
of [23]. Fluctuations in concentration ((><!>) and shear rate 
(5j) evolve via the diffusion and Navier-Stokes equations, 
in which shear-induced migration and the <!> dependence 
of the shear stress must be included [23] . For small fluc- 
tuations along z we have (to linear order in Sj): 

d t 5$ = -V • J ~ M (n* <9 2 2 (5$ + d 2 Jj) , (1) 
d t Sj = p- x d 2 z a ~ p- 1 9f<5$ + 8*5j) . (2) 

Here we have introduced the shorthand = dU/djl^,, 
likewise n$ , o~j and 0$ ; M is a collective mobility and 
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FIG. 3: Line: critical flow rate Pec^) from Eq. J7J) with 
[p,r,m = n, A, B] = [3,4,0.4,25,0.005]. Colors mark the 
value of F. Symbols: experimental values for Pe c (averaged 
over z g ). Arrow: possible evolution of an unstable state. 



p the density. The migration current J arises from par- 
ticle pressure gradients 9 Z II($,7) due to variations in 



both <E> and 7 



|35| . The terms involving II j and 
(T$ in Eqs. (|1I2I) cause respectively particle migration to- 
wards regions of low shear rate, and accelerated shear 
in regions of low concentration; together, these amplify 
fluctuations. That is, a fluctuation towards higher <£> in 
some region creates a lower shear rate there. This pro- 
motes inward particle migration, giving a positive feed- 
back effect. This tendency is counteracted by the re- 
maining terms which describe stable diffusive spread- 
ing of both particles and momentum (or equivalently 
shear rate). Rewriting Eqs. (|U2j) as dt^i = Lijdl^j 
with — (5Q,8j), we see that instability sets in when 
det Lij = M(n$er.y — n^crg>)/p becomes negative, or 



equivalently, as first derived in 



when 



F = 



_ j.x 7 



> 1. 



(3) 



To evaluate F, we first write the HB form for a in terms 
of the Peclet number Pe = 77-5 : 

The first term is a y and n ~ 0.4 — 0.5, as before. Typical 
values for A from our experiments are A = 10 — 20. The 
particle pressure IT for HS colloids has a similar form: 

n $ 



n 



(1-$) 



[l + .g($)Pe m ], ff ($) = B(l - $) 



(5) 



Here ITo = 2.1750 m fcs7 1 '/ 'ira 3 ; the first term then ap- 
proximates the osmotic pressure at rest [36| , whereas the 
second term is the contribution due to shear. For the 
latter, a Pe 2 dependence was found in (3(| but this 
is restricted to a linearly viscous regime (cr oc Pe) at very 
small Pe, in which Eq. [3] yields F oc Pe 2 resulting in sta- 
ble flow. However, this regime is hard to access, and in- 
deed completely absent whenever a y is nonzero ($ > $ g ), 
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FIG. 4: Velocity profiles v(z)/v(z g ), calculated from 7 oc 
[{a/ay ($(«))) - l] 2 and a y _oc (1 - $)~ 3 , with $(z) - $ = 
\6$\[{z/z g ) - (1/2)] with $ = 0.96, |(5$|/$ = 0.002, for 
different <j/<t h ($). (b) Effective flow curves a versus 7 = 
[v(z g , cr, $(z))]/z g for various $ and \5$\. 



where it is replaced by a non-Newtonian regime in which 
both the stress and the particle pressure increase sublin- 
early with Pe. From simulations for $ < <J> 3 [25|, |37j , we 
extract to = 0.4 — 0.5, (and B ~ 0.003, r ~ 3), so that 
to — 71, as is also observed in simulations of 2D foams 
[38j . For $ > $ 5 , and for glassy flow in general, we ex- 
pect these two exponents to remain similar (e.q^, within 
mode coupling theory, m = n seems probable [20j). 

Using Eqs. (|4l5p . we obtain a limiting value for F at 
large Pe as F^, — m(p — n)/nr 39]. The flow is thus 
stable in this regime provided r > p — n; to explain our 
experiments we require r > 2.5. More interesting is the 
result for small Pe, where we obtain 



F -> Pe" 



-n mpg{<S>)<S> 



mpB<& 



ns($) nA(l - §) r + r - 



(6) 



Here Fq(<&) is a quasi-plateau value maintained while 
Pe m_ ™ ~ 1 (and a true limiting value if to = n). It 
follows that homogeneous flow at low Pe is unstable for 
concentrations $ > <1> C , where Fo($ c ) = 1- We argue 
that it is this SCC-induced instability that creates the 
small variations <5$(z) that were assumed in Fig. HJa), 
and which account for the experimentally observed lo- 
calized flow at small 7. Because, as in [23j, our analysis 
is limited to linear stability, we do not have a clear idea 
of the mechanism limiting the growth of 5$. However, 
the extreme nonlinear dependence of both a and n on $ 
as $ — » 1 makes it plausible that <5$ remains small. 

The critical flow rate Pe c , below which the instability 
sets in, follows from the condition f 1 (Pe c ,<I ) ) = 1. For 
the case m — n this simplifies to j39j : 



Pe c ($) 



P 



rs($)(l-^°°) 



1 - 



1 



F ($) 



l/n 



(7) 



Figure [3] shows Pe c ($) for specific parameter values. En- 
tering the unstable regime $ > Q c , Pe c increases and 
grows ~ s($) -1 / n for $ — > 1. Our model is able to 
give a quantitative account of the data; the fit [4(J gives 



4 



$ c = 0.89, close to, but above, $ g . (Note that Eq. © 
does not rule out $ c < $ 9 in principle.) 

In contrast to most shear banding scenarios the 
present system seems unable to achieve global stability 
by separating into distinct bands. This is illustrated by 
the arrows in Fig. [3] an initially unstable state can form 
a locally depleted region, $(z) < $, that is stabilized 
when Pe($(z)) > Pe c , but the remaining concentrated 
region is even more unstable than before. This suggests 
that the banded flow should have residual temporal fluc- 
tuations. We have indeed observed this (Fig.fTJd)); after 
startup shear of 7 = 0.01 s _1 , a weakly nonlinear profile 
develops a central region where j(z) is strongly reduced, 
which then speeds up by expanding the lower band, and 
subsequently reverts to a larger j(z). We have also ob- 
served (in a planar shear cell [3l[) that a fluidized band 
can swap from z ~ to z ~ z g over sufficient time. 
Leaving a detailed study of these effects for future work, 
we note that the shear banding we observe is permanent 
(if unsteady) not transient; moreover, our stress is time- 
independent, in contrast to the results in (T7j . 

Our shear-banding mechanism arises from the con- 
centration dependent nonlinear rheology of glasses along 
with the nonlinear process of 'Brownian dilation' whereby 
flow increases the particle pressure n($,Pe). For larger 



Pe than those studied here, hydrodynamic effects lead 
to a much stronger, linear increase of II with flow rate 
[26| . ultimately leading to shear thickening in these sus- 
pensions. Hence, the localization we observe might be 
interpretable as a precursor to shear thickening. 

In conclusion, we have shown that HS colloidal glasses 
exhibit a new type of shear banding, well described by a 
model in which SCC leads to unstable flow near yielding. 
In this scenario, very small concentration variations can 
sustain large variations in flow rate. Our results may 
also be relevant for flow in other glassy materials, such as 
foams where 'dilatancy' has recently been observed [4~L| . 
It may also shed light on shear bands and dilatancy in 
metallic glasses [5|. We mention in particular simulations 
of a model glass [18] , which showed both fluctuating shear 
bands and small but finite concentration gradients. 
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